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We present numerical results for the 2-flavour Schwinger model with dy- 
namical chiral lattice fermions. We insert an approximately chiral hy- 
percube Dirac operator into the overlap formula to construct the overlap 
hypercube operator. This is an exact solution to the Ginsparg- Wilson re- 
lation, with an excellent level of locality and scaling. Due to its similarity 
with the hypercubic kernel, a low polynomial in this kernel provides a nu- 
merically efficient Hybrid Monte Carlo force. We measure the microscopic 
Dirac spectrum and discuss the corresponding scale-invariant parameter, 
which takes a surprising form. This is an interesting case, since Random 
Matrix Theory is unexplored for this setting, where the chiral condensate 
E vanishes in the chiral limit. We also measure E and the "pion" mass, 
in distinct topological sectors. In this context we discuss and probe the 
topological summation of observables by various methods, as well as the 
evaluation of the topological susceptibility. The feasibility of this summa- 
tion is essential for the prospects of dynamical overlap fermions in QCD. 
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1 The Schwinger model 

The Schwinger model describes QED in two dimensions, i. e. 2-component 
Dirac fermions interacting through a U{1) gauge field [1]. In a Euclidean 
plane, the Lagrangian reads 

It is a popular toy model for QCD; in particular it shares the property of 
confinement [I1I2]- On the other hand there are fundamental differences, 
such as the super-renormalisability of the Schwinger model, and a non- 
running gauge coupling g. 

A further qualitative difference — which is of particular interest in this 
work — is the spontaneous chiral symmetry breaking in QCD with massless 
quarks. In d = 2 this effect can be mimicked to some extent for instance 
by the Gross-Neveu model, where a discrete chiral symmetry breaks spon- 
taneously. However, in the Schwinger model with fermion mass m = the 
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chiral symmetry is continuous, and therefore it cannot undergo sponta- 
neous symmetry breaking (SSB) due to the Mermin- Wagner Theorem [5]. 
Nevertheless the chiral condensate E = —{iptp) , which acts as the order 
parameter for chiral symmetry breaking, takes a non-vanishing value in the 
1-flavour case, because of the explicit symmetry breaking due to the axial 
anomaly pQ. This leads to S(m = 0) = [e' /2tx'^/'^) g c:^ 0.16(7 (where 7 is 
Euler's constant). 

For Nf > 2, however, the massless limit has an unbroken chiral symme- 
try. For Nf degenerate fermion flavours of mass m the chiral condensate 
behaves as |1] 



EMoc(— j ^ * = (1-2) 

where /3 = 1/g'^, and 6 is the critical exponent. 

In our study we consider Nf = 2. Here there are analytical evaluations 
of the proportionality constant for the case of light fermions, m <^ g, based 
on bosonisation and low energy assumptions, 

, , fmy/^ ^ / 0.372... Ref. [5] ,^ 

S(m) = const, (-j , const. = | ^ ^^g _ _ _ . (1.3) 

Under the same assumptions the mass of the iso-triplet ( "pion" ) ^ and of 
the iso-singlet ("■?] particle") [7] are predicted as 



= 2.008 . . . {m^gY^^ , = \j . (1.4) 

The relation M^r oc m^/^ replaces the Gell-Mann-Oakes-Renner relation 
(M^ oc y^) of QCD [8J. 

For models with a finite condensate E(m = 0), its value can be deter- 
mined from the Banks-Casher plateau of the Dirac eigenvalue density at 
zero |9] (or near zero in a finite volume). Moreover, chiral Random Matrix 
Theory (RMT) has elaborated subtle techniques to predict a wiggle struc- 
ture on this plateau, which allows for a refined determination of S from 
the densities of the low-lying Dirac eigenvalues [TOl|Tl]. This method has 
been tested successfully in the e-regime of QCD with quenched [T2]4T^ 
and with dynamical [Ml [17] quarks, and also in the 1-flavour Schwinger 
model [T8t[T9]. The latter studies were based on configurations, which were 



^Strictly speaking E diverges logarithmically for increasing volume in the quenched 
approximation [15]. Still quenched QCD in boxes of length >1.2 fm yields sensible 
results for the chiral condensate. 



3 



generated quenched and later re- weighted with the fermion determinant; 
we denote this method as "quenched re- weighted" H 

However, the estabhshed RMT techniques are not apphcable in the 
chiral hmit of the 2-flavour Schwinger model; RMT for situations with 
E(m = 0) = awaits to be worked out. Nevertheless we confirm the usual 
RMT prediction for the unfolded level spacing distribution in a unitary 
ensemble. On the other hand, the microscopic spectrum does not exhibit 
a Banks-Casher plateau. Instead we observe to high precision the scale- 
invariance of the product XV^/^^ where A is a low-lying Dirac eigenvalue 
in the volume V . This result remains to be understood from the RMT 
perspective, since it cannot be explained simply with the critical expo- 
nent 5 = 3 given in eq. fll.2p . We also discuss the densities of the Dirac 
eigenvalues in the bulk and their scaling behaviour. 

Next we confront the measurement of S(m), based on the full Dirac 
spectrum, with eq. (11. 3p . This requires an (approximate) summation of 
the values in all topological sectors, guided by the measurements in a few 
sectors. We probe several methods for this purpose and apply them to 
E and to the "pion" mass given in eq. (II. 4p . These approaches also in- 
volve a determination of the topological susceptibility. The requirement of 
a topological summation is generic for simulations with dynamical chiral 
fermions, because the Monte Carlo histories tend to perform only very few 
topological transitions. Also other lattice fermion formulations, such as the 
Wilson fermion and variants thereof, will run into the same problem when 
they represent light fermions on very fine lattices |22] (say a < 0.05 fm in 
QCD). Therefore the applicability of these techniques is relevant, particu- 
larly in view of QCD simulations with light quarks close to the continuum 
limit. 

Section 2 describes our lattice formulation of the Schwinger model and 
discusses some of its properties, in particular locality and scaling. Section 
3 presents our version of the Hybrid Monte Carlo algorithm which we used 
in this study. We discuss its properties regarding conceptual conditions, 
and practical aspects of its performance. Section 4 deals with the Dirac 
spectrum, the construction of a scale-invariant variable and the link to 
RMT. Section 5 discusses the summation of S and over the topological 
sectors, along with the evaluation of the topological susceptibility. Sec- 
tion 6 is devoted to our conclusions. Progress reports of this project have 
appeared in several proceeding contributions [23] . 

^That method was pioneered in the Schwinger model in Ref. [20 . It worked success- 
fully in some cases, but it runs into trouble when the fermion determinant fluctuates 
strongly, as it happens for very light fermions. Hence the study presented here requires 
the simulation of truly dynamical fermions, as it was first attempted in Ref. |21) . 
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2 Lattice formulation with overlap Hyper- 
cube Fermions 



We consider the lattice formulation of the Schwinger model with compact 
link variables U^^fj, G f/(l), and with the plaquette gauge action. Remark- 
ably, for the pure gauge theory this is indeed a perfect lattice action 
For the fermions we employ an overlap hypercube fermion (overlap-HF) 
Dirac operator, which is an exact solution to the Ginsparg- Wilson Rela- 
tion (GWR). 

The GWR is a criterion for a lattice modified, exact chiral symme- 
try |25], which first emerged from the study of perfect actions for lat- 
tice fermions [2 ^ [261428] . Independently, chiral lattice fermions were con- 
structed in the Domain Wall Fermion formulation [29] , which separates the 
zero modes of opposite chirality in an extra "dimension" H Integrating out 
this extra direction leads to the overlap formula [SOj, which provides yet 
another way to represent a chiral vector theory on the lattice [SI]. The 
lattice Dirac operator for Domain Wall Fermions (in the limit of an infinite 
wall separation) and for overlap fermions turned out to be solutions to the 
GWR as well [3T]. 

Its importance as a general chirality criterion was first pointed out in 
Refs. [2B] , which showed in addition that classically perfect fermion actions 
obey this criterion as well. Since those formulations involve couplings over 
an infinite range (in > 1), a truncation is needed, which marginally dis- 
torts the perfect symmetry and scaling properties. For the free, optimally 
local, perfect fermion [2l] the truncation to a unit hypercube by means of 
periodic boundary conditions over three lattice spacings preserves excellent 
scaling [52] and chirahty [33]. It leads to a lattice Dirac operator of the 
form 

DnF{x,x + r) = p^{r)-f^ + \{r) , (2.1) 

i. e. a vector term plus a scalar term, as in the case of the Wilson fermion, 
but with an extended structure {x and x + r are lattice sites). 

In d = 2 these terms include only couplings to nearest neighbour lattice 
sites and across the plaquette diagonals. We are using here the version 
denoted as CO-HF (Chirally Optimised Hypercube Fermion) in Ref. |34j . 
which is optimal for our algorithm to be described in Section 3. For con- 
venience we display in Tabled] the couplings in the notation of eq. fl2.ip . 

We gauge Dhf by multiplying the compact link variables along the 
shortest lattice paths connecting x and y = x + r ; for the diagonal the two 
shortest paths are averaged [51]. Thus we arrive at the operator Dup^xyiU), 

■^An extra direction is introduced, which appears as a dimension for the free fermion, 
but which does not carry gauge fields. 
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r 


Plir) 


A(r) 


(0,0) 





1.49090692 


(1,0) 


0.30583220 


-0.24771369 


(1,1) 


0.09708390 


-0.12501304 



Table 1: The coupling constants of the Chirally Optimised Hypercuhe 
Fermion (CO-HF) 1341- Note that pi(r) is even in ri and odd in r2 (and 
vice versa for P2{j')), while A(r) is even in both components ofr. 

which characterises the interacting Hypercube Fermion (HF). 
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Figure 1: The spectra of D^p and of Dovuf (with and without mass) in 
the complex plane, for a typical configuration generated at P = 5 on a 
16 X 16 lattice with m = 0.01 (on the left) and at m = 0.24 (on the right). 
The similarity of D^p to D^^^^ and to Do^i{p{m) shows that the HF is 
approximately chiral, and useful for an efficient computation of the Hybrid 
Monte Carlo force (see Section 3). 

Since the operator D^p is "73-Hermitian" , Dl^p = 73£'hf73, the exact 
chirahty (which got lost in the truncation) can be restored by inserting Dhf 
into the overlap formula [31]. This yields the overlap-HF operator j33|[33] 

^ovHF(m) = (1 - ^) ^ivHF + ^ , 

DZp = 1+73-^, i^HF = 73(^HF - 1) . (2.2) 
V^HF 

i^HF is Hermitiai]0 and i^o^HF fulfils the GWR in its simplest form, 

WZf,!^} = D^F^sD^p . (2.3) 

^The constant 1 in the formulation of H^p is fine since we deal with smooth gauge 
configurations. For stronger gauge couplings one would prefer to increase this constant. 
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In practice we evaluate this overlap operator by means of rational Zolotarev 
polynomials, as suggested in Ref. after projecting out the lowest two 
modes of Dl^-pDuF, which are treated separately. 

Compared to H. Neuberger's standard overlap operator [S], we 
replace the Wilson kernel I^w by Drf [33]. Since the latter is an approxi- 
mate solution to the GWR already, the transition D^p -DqvHF is only a 
modest chiral correction, 

Do.HF ~ DnF , (2.4) 

in contrast to the drastic transition — )■ D^. This property is illustrated 
in Figure [H which compares the spectra of Duf and -DqvHF for typical 
gauge configurations generated at m = 0.01 and at m = 0.24, both at 
/3 = 5 on a 16 X 16 lattice. The similarity of D^p to -Dq^hf useful for 
the computation of the overlap operator and for its favourable properties 
in addition to chirality (see below), while the similarity to DovHf('^) is 
helpful for our algorithm to be discussed in Section 3. 

Due to its perfect action background, Dhf is also endowed with a good 
scaling behaviour and approximate rotation symmetry, which are inherited 
by -DqvHF thanks to relation f l2.4p . That relation further provides a high 
level of locality for -Dq^hf ' since it deviates only little from the ultralo- 
cal operator Dhf • All these properties have been tested and confirmed 
extensively in the quenched re- weighted study of Ref. [M] . 

Let us reconsider here the level of locality, which is a key criterion in 
the comparison of different chiral lattice fermion formulations. We test it 
in the usual way , by applying -Dq^hf ^ ^'^i^ source rj and measuring 
the decay of the function 

/(r) = T {DZF,.yiU) Vy I E l^/^-^^l ='}^ Vy = Syo(l) ■ (2.5) 

We first consider the free fermion and demonstrate that this decay is clearly 
faster for the overlap-HF operator -DqvHF than for the Neuberger operator 
Dtss, see Figure |2]on top. The plot below shows that the decay is still expo- 
nential for the configurations that we generated with dynamical fermions 
at /3 = 5, which confirms the locality of our Dirac operator. This assures 
that our lattice fermion formulation is conceptually on safe grounds. In 
the range that we studied, the mass has practically no influence on this 
decay rate. We observe that -DqvHF has a higher degree of locality than 
Z^N, since -DqvHF at /3 = 5 is still clearly more local than even the free : 
the decay for the free -DqvHF, f(j) oc exp(— 1.5r), is reduced just slightly to 
exp(— 1.45r) by the gauge interaction, whereas only decays as exp(— r) 
even in the absence of gauge fields. 
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This virtue also holds for the overlap-HF in quenched QCD [HIEZ]: at 
/3 = 6 the exponent is increased by almost a factor of 2 compared to Dn, 
and the locality of overlap-HF is manifest down to /3 = 5.6. This enables 
the formulation of chiral fermions on coarser lattices than the use of Dn, 
which is of importance in view of QCD simulations at finite temperature. 
In that case, simulations are performed with a very small number Nf of 
lattice sites in the Euclidean time direction. Its extension is extremely 
expensive; the computational effort grows oc N^'^ [38]. The application of 
the — and in future also of -DqvHF — is therefore most promising in 
finite temperature QCD [39] . 

The scaling behaviour was found to be excellent for both, the HP and 
the overlap-HF, by considering dispersion relations in the free case and in 
the 2-flavour Schwinger model with quenched re-weighted configurations, 
which were generated at /3 = 6 [31]. The HP and the overlap-HF have 
an even better scaling behaviour than the (truncated) classically perfect 
action, which was constructed and tested for the Schwinger model in Ref. 
[20] (although that concept was actually designed for asymptotically free 
theories [10])- Another quenched re- weighted scaling test was added in 
Ref. [H]. 

Throughout this work we fix /3 = 5 and study the effects of varying 
the lattice size and the fermion mass. So we do not investigate explicitly 
the continuum extrapolation, since the scaling artifacts due to the finite 
lattice spacing turned out to be very small. This is illustrated by the 
dispersion relations of the "meson" masses shown in Figure [31 they follow 
the continuum behaviour up to quite large momenta, much further than the 
Wilson fermion or the Neuberger fermion [Mj. Moreover scaling artifacts 
are expected to be negligible also based on the large plaquette values near 
0.9, see Table ID Hence the configurations are smooth, which corresponds 
to a fine lattice. 

On the other hand, the issue of finite size effects is relevant here, and 
we will address it extensively in Sections 4 and 5. Figure H] shows the 
correlation length ^ = 1/M^ as a function of the fermion mass, as expected 
in infinite volume according to eq. f ll.4p . It reveals that significant finite 
size effects may occur for our smallest fermion masses and volumes. These 
effects can be very useful to investigate the distinction between topological 
sectors. In QCD they have been used to determine some of the Low Energy 
Constants by means of simulations in — or close to — the e-regime [T2] - 
[n ] [TB ] [T7 ]li2] - |il] and the 5- regime In our study the finite size effects 
are useful since they provide a suitable laboratory to probe methods of 
summing up observables measured separately in a few topological sectors. 
Section 5 presents pilot studies of such procedures, which might become 
relevant in lattice QCD. 
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Figure 2: The locality of the overlap Dirac operators, tested by the decay of 
the function l{2.5\) . against the taxi driver distance r = |ri| + |r2| . On top 
we compare our overlap-HF operator D^^p to Neuberger's standard overlap 
operator D-^ (with a Wilson kernel) in the free case. ^4^ r = 3 (r = 7), /(r) 
is already suppressed by more than one (two) order (s) of magnitude for 
-^ovHF- P^^^ below shows the exponential decay of (/(r)) based on our 
overlap-HF simulations at (3 = 5. The gauge interaction reduces the decay 
rate just marginally, with hardly any dependence on the fermion mass m. 
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Figure 3: The dispersion relations for the "pion" and the '^-particle" 
(isospin triplet and singlet), measured at fermion mass m = 0.01 at L = 20 
and L = 32, in different topological sectors. Irrespective of the small dif- 
ferences, they follow in all cases very closely the theoretical curve in the 
continuum, given by eq. Up to momentum p ^ n/2 lattice artifacts 

are tiny, which confirms an excellent scaling behaviour. This scaling qual- 
ity is similar to the overlap-HF in our previous quenched re-weighted study 
at P = Q, but in that case the Wilson fermion and the Neuberger fermion 
show sizable scaling artifacts setting in at p ^ 0.9 JS^ . 
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Figure 4: The correlation length ^ = (in infinite volume) as a func- 

tion of the degenerate fermion mass m, according to eq. For instance 

for our lightest fermion mass, m = 0.01, it amounts to C, — 14. 

3 Hybrid Monte Carlo with a preconditioned 
force 

3.1 Algorithm 

In order to simulate overlap-HFs dynamically, the standard Hybrid Monte 
Carlo (HMC) algorithni^ would use the fermionic force term 

7n-l fn-^ ^QovHF I ^QovHF ^-1 / /q 1 \ 

y ^ovHF [yovHF r VovHF J VovHF V j W-^) 

where QovHF = 75-DovHF is the Hermitian overlap-HF operator, and A^^fj, 
are the non-compact gauge link variables. However, this force term is 
computationally expensive. In particular in view of prospects for QCD we 
are going to explore a hopefully efficient alternative. In addition, the force 
(13.11) is conceptually problematic due to the discontinuous sign function 
Hrf I a/ -f^HF QovHF, see eq. (12. 2p . Under a continuous deformation of 
the configurations, the denominator vanishes at the transition between 
topological sectors. Here we refer to the proper definition of the topological 
charge by means of the fermionic index ^28j 

1/ = n_ — n+ , (3.2) 

^The original work on the HMC algorithm is Ref. [46^; pedagogical descriptions can 
be found for instance in Refs. [4 7) . 
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where ?t,+ is the number of zero modes of -DqvHF with positive (neg- 
ative) chirahty|§ Regarding the spectrum of Hup, a topological transition 
means that an eigenvalue crosses 0, so its map by the overlap formula f l2.2p 
flips between m and 2, i.e. it either appears as a zero mode of -DqvHF ' 
it is sent to the cutoff scale. Indeed the QCD simulations that have been 
performed with this HMC force hardly ever achieve topological transitions, 
see in particular Refs. [HlllH]; this issue is discussed in detail in Ref. [i9] . 
We repeat that the same problem occurs also with other lattice Dirac oper- 
ators for light dynamical fermions when the lattice becomes very fine [22] : 
the HMC history tends to get trapped in one topological sector. 

Here we report on HMC simulations, which are again facilitated by the 
powerful property (12. 4p . Our algorithmic concept follows the HMC version, 
which was applied to improved staggered fermions of the HF-type in Ref. 
[50] . It used a sophisticated Dirac operator in the accept/reject step, but 
a simplified formulation for the force, which is quick to evaluate. In order 
to simulate the dynamical overlap-HF we render the force term continuous 
and computationally cheap by inserting only approximate overlap operators 
in the term (13. ip . To be explicit, we apply an overlap-HF to a low precision 
e' in the external factors QovHF' we use Huf instead of QovHF in the 
derivatives (although this could be extended to a polynomial as well)Jll 



For the integration we applied the Sexton- Weingarten scheme [5T] with a 
partial (St)^ correction (where St is the step size). The time scales for the 
fermionic vs. gauge force had the ratio of 1 : 5, but we did not observe much 
sensitivity to this ratio. 

The Metropolis accept/reject step is performed with the high precision 
overlap operator -DqvHf.e- Hence the deviations in the force are corrected, 
and the point to worry about is just the acceptance rate. The complete 
simplification, which reduces QovHF,e' to 75DHF, was considered in Ref. [52] . 
which found a decreasing acceptance rate for increasing volume; that study 
was based on the Scaling Optimised Hypercube Fermion (SO-HF) of Ref. 
[31]. However, in this respect it turns out to be highly profitable — and 
still cheap — to correct the external factors to a modest precision. We 

^We fix the sign of tfie index sucli tliat it niatclies tlie continuum gauge formulation 
of the topological charge, j (Px ei2Fi2- However, in all considerations of this work only 
matters. 

^The simplified force that we are using is not only computationally cheaper, but it 
is also expected to be helpful to achieve at least a few topological transitions. Hence it 
might not even be favourable to extend -ffuF in eq. p.3p to a polynomial, which would 
move the force formulation closer to the dangerous sign function. 




(3.3) 
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m 


number of configurations 


topological 




iy = 


|z/| = 1 


|z/| = 2 


total 


transitions 


0.01 


2428 


307 




2735 


7 


0.03 


1070 


508 




1578 


2 


0.06 


741 


660 




1401 


7 


0.09 


919 


587 


1 


1507 


7 


0.12 


664 


501 


248 


1413 


8 


0.18 


791 


563 


50 


1404 


15 


0.24 


576 


978 


56 


1637 


17 



Table 2: Our statistics of configurations at L = IQ and seven fermion 
masses m. The HMC trajectory lengths i are given in Table \^ In all 
cases they consist of 20 integration steps (6t = i/20). The measurements 
are separated by at least 200 trajectories. For m = 0.01 this separation 
was enlarged to 600 trajectories for a better decorrelation. As a generic 
property of dynamical overlap fermion simulations, topological transitions 
are rare, as we see in particular for our lightest fermion masses. 

chose the algorithmic parameters for the (absolute) precisions as 

e' = 0.005 (force term) , e = lO'^^ (Metropohs step) , (3.4) 

which increases the acceptance rate by an order of magnitude compared to 
the simple use of H^p throughout the force term. The force we obtain in 
this way is not based on Hamiltonian dynamics, but the way it deviates 
from it (by proceeding from 75i^HF to QovHF.e') does maintain the property 
of area conservation in phase space. 

3.2 Statistics 

With this algorithm, we performed production runs at /3 = 5 on L x L 
lattices of the sizes L = 16, 20, 24, 28 and 32. On the 16 x 16 lattice we 
collected statistics at seven fermion masses in the range m = 0.01 . . . 0.24 in 
distinct topological sectors, as displayed in Table |2i At our lightest mass, 
m = 0.01, we performed additional simulations on larger lattices of size 
L = 20 ... 32, plus runs at m = 0.06, L = 32, see Table O 

Since the force term (13. 3p tends to push the trajectory a bit off the 
hyper-surface of constant energy, we kept the trajectory length i (between 
the Metropolis steps) short. At L = 16, m < 0.18, we chose i = 1/8, which 
is divided into 20 integration steps {6t = 0.00625). This was a compromise 
in view of precision and dynamics between the trajectory end-points. On 
the larger lattices the trajectory length i was further reduced, see Table 
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L 


m 


z/ = 


number 
|z/| = 1 


of configi 

\v\ = 2 


nations 
W\ = 3 


total 


20 


0.01 


435 


304 






739 


24 


0.01 






278 


273 


551 


28 


0.01 




240 




180 


420 


32 


0.01 


138 


98 


82 




318 


32 


0.06 


91 


293 






384 



Table 3: Our statistics for the lattice sizes L = 20 ... 32 at fermion masses 
m = 0.01 and 0.06. The trajectory length i was reduced for increasing L, 
see Table^ while the integration step was always fixed as 6t = i/20. After 
thermalisation, the configurations are separated by at least 200 trajectories. 
We show our statistics in the topological sectors |z/| = . . . 3. In particular 
at m = 0.01 topological transitions were very rare, so we captured various 
sectors by a multitude of "hot starts". 

ml 

The configurations used for the measurements were separated by at least 
200 trajectories. Still the autocorrelation with respect to the observables 
in Sections 4 and 5 is not always negligible, see Table [51 In particular 
some problems show up for L > 20 and higher topological charges, which 
suggests that an application of this algorithmic approach to QCD might 
require further refinements. Here autocorrelations were taken into account 
by a jackknife analysis of the measured data. That method was used for 
the error calculations throughout this work; we probed a variety of bin 
sizes and took the maximal error in each case. 

3.3 Reversibility, acceptance rate and computational 
effort 

Reversibility — to a good precision — is a requirement of the HMC algo- 
rithm. The possible danger for this crucial property could be an instability 
in the molecular dynamics trajectory due to directions with positive Lya- 
punov exponents, such that certain deviations from the exact trajectory 
are amplified exponentially [53] . 

To test if we are confronted with this problem, we moved forth and back 
with a variable number of steps, and measured the (absolute) modification 
of the gauge action, lAS^I- As examples. Figure O (on top) shows our 
results for the reversibility precision at L = 16, 6t = 0.00625, for the 

^At a few points in the HMC histories, where the algorithm run into convergence 
problems, we temporarily reduced i below the trajectory length given in Tabled always 
maintaining the dissection into St — £/20. 
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masses m = 0.03, 0.12 and 0.24. The precision is satisfactory in all cases. 
It still improves significantly for increasing mass, as we also observe for 
6t = 0.005, see Figure |5] (below). Our results do not indicate the presence 
of any positive Lyapunov exponent. 




15 20 25 30 ;i5 

number of steps of size St = 0.00625 




number of steps of size St = 0.005 



Figure 5: The reversibility precision with respect to the gauge action Sq for 
a variable number of steps of length 6t = 0.00625 (on top) and 6t = 0.005 
(below), both on a 16 x 16 lattice at 13 = 5. We show results for very different 
masses. There is no indication for any positive Lyapunov exponent. The 
precision improves as we increase m, but it is satisfactory in all cases. 

Being confident that our algorithm is safe, we now proceed to the ques- 
tion of its efficiency. The plots in Figure [6] show the acceptance rates in 
the Metropolis step@ In some sectors they were somewhat modest for the 
parameters chosen here, which is related to the aforementioned cases of 
rather long autocorrelation times. 

^We evaluated the acceptance rate by considering the acceptance probabihty in each 
MetropoUs step, regardless of the actual accept/reject decision. This is statistically 
more conclusive than just counting the acceptance ratio. 
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Figure [7] displays the number of Conjugate Gradient iterations that 
was required per trajectory, specifically in the evaluation of -DqvHF (upper 
plots) and in total (lower plots). Table H] summarises the acceptance rates, 
as well as the number of Conjugate Gradient steps in the evaluation of the 
overlap operator and in total. We add the plaquette value to characterise 
the smoothness of the configurations; this serves as a point of orientation 
for comparison with other models in lattice gauge theory. 
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Figure 6: The Metropolis acceptance rate on the L = 16 lattice at seven 
different masses (on the left), and at m = 0.01 on five lattice sizes (on the 
right). Note that the trajectory length varies, as specified in Table\^ 



As usual, the simulation becomes faster when we proceed from light 
to moderate fermion mass. However, the number of Conjugate Gradient 
iterations rises only mildly as we approach the chiral limit, even down to 
very light masses. The reason is that finite size effects prevent the non-zero 
eigenvalues from becoming really tiny (this virtue is obviously reduced if 
the volume increases). The low- lying non-zero eigenvalues will be discussed 
in detail in the next section. 



4 The Dirac spectrum 



In this section we discuss our results for the eigenvalues of the Dirac opera- 
tor -Dqvhf ^l- (l2-2p . after stereographic projection (a Mobius transform) 
onto the half-line IR 



A,; 



1 - Xi/2 



(4.1) 



4.1 Unfolded level spacing distribution 

We first consider the full spectra and the resulting unfolded level spacing 
distribution. This distribution is obtained as follows: one first numerates 
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L 


m 


i 


acceptance 


number of CG iterations 


plaquette 








rate 


in i^ovHF 


total 


value 


16 


0.01 


0.125 


0.439(8) 


97.2(1) 


4529(4) 


0.8971(2) 


16 


0.03 


0.125 


0.295(15) 


94.2(1) 


4419(5) 


0.8965(4) 


16 


0.06 


0.125 


0.306(17) 


93.9(1) 


4405(4) 


0.8974(4) 


16 


0.09 


0.125 


0.355(9) 


91.4(1) 


4312(2) 


0.8963(2) 


16 


0.12 


0.125 


0.386(12) 


89.6(1) 


4244(3) 


0.8961(3) 


16 


0.18 


0.125 


0.450(9) 


82.7(1) 


3966(2) 


0.8965(2) 


16 


0.24 


0.0625 


0.706(13) 


76.5(1) 


3708(3) 


0.8947(4) 


20 


0.01 


0.0625 


0.504(13) 


128.7(2) 


5961(6) 


0.8972(3) 


24 


0.01 


0.05 


0.352(9) 


149.9(6) 


7129(27) 


0.8969(5) 


28 


0.01 


0.04 


0.242(13) 


191.9(2) 


9152(17) 


0.8972(5) 


32 


0.01 


0.03 


0.342(16) 


224.9(5) 


10432(18) 


0.8968(5) 


32 


0.06 


0.03 


0.619(10) 


199.0(2) 


9546(8) 


0.8971(3) 



Table 4: The characteristic indicators for the performance of the precondi- 
tioned HMC algorithm: first we give the acceptance rate in the Metropolis 
step at the end of each trajectory of length I. The accept /reject step uses 
-DovHF to machine precision (see eq. Ii3.4^ ). We add the number of Con- 
jugate Gradient iterations needed to evaluate the operator -DqvHF d'nd in 
total. Finally we quantify the smoothness of the configurations by the mean 
plaquette value. 
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L 


m 


z/ 


plaquette 


Dirac eigenvalue Ai 


chiral condensate E 


16 


0.01 





0.6 


0.7 


1.1 


16 


0.01 


1 


0.7 


0.5 


0.8 


16 


0.03 





0.6 


1.2 


1.2 


16 


0.03 


1 


0.5 


1.0 


1.1 


16 


0.06 





0.6 


1.2 


1.2 


16 


0.06 


1 


0.7 


0.9 


0.9 


16 


0.09 





< 0.5 


1.2 


1.2 


16 


0.09 


1 


< 0.5 


0.9 


1.0 


16 


0.12 





< 0.5 


1.5 


1.3 


16 


0.12 


1 


0.5 


0.9 


0.9 


16 


0.12 


2 


0.5 


0.8 


0.9 


16 


0.18 





< 0.5 


1.1 


0.9 


16 


0.18 


1 


0.6 


0.7 


0.7 


16 


0.18 


2 


< 0.5 


0.6 


0.6 


16 


0.24 





< 0.5 


1.1 


1.1 


16 


0.24 


1 


< 0.5 


0.9 


0.7 


16 


0.21 


2 


< 0.5 


0.5 


< 0.5 


20 


0.01 





< 0.5 


1.3 


0.9 


20 


0.01 


1 


< 0.5 


1.2 


1.1 


24 


0.01 


2 


< 0.5 


2.5 


2.5 


24 


0.01 


3 


0.6 


2.8 


2.4 


28 


0.01 


1 


0.5 


2.7 


1.7 


28 


0.01 


3 


0.5 


4.5 


5.0 


32 


0.01 





< 0.5 


1.3 


1.2 


32 


0.01 


1 


< 0.5 


1.8 


1.6 


32 


0.01 


2 


1.3 


1.7 


1.9 


32 


0.06 





< 0.5 


0.9 


0.8 


32 


0.06 


1 


< 0.5 


1.9 


2.7 



Tabic 5: The integrated autocorrelation times Tint = | + Sr>i ^i'^) (where 
C{t) is the autocorrelation function) over a total trajectory length 25, for 
the mean plaquette value, the leading non-zero Dirac eigenvalue Ai ( relevant 
in Section 4), and the chiral condensate E (relevant in Section 5). 
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Figure 7: The number of Conjugate Gradient iterations per trajectory in 
DovUF (upper plots), and including all operations (lower plots). We show 
results at L = IQ and various masses (plots on the left), and at m = 0.01 
on various lattice sizes (plots on the right). 



the eigenvalues of single configurations in ascending order, Aj (here we omit 
the eigenvalues with negative imaginary parts before the mapping (14. ip ). 
Next we put all eigenvalues in a set of configurations together and numerate 
them again. Now we consider pairs of eigenvalues Aj, Aj+i (of the same 
configuration), and count by how many ordinal numbers they differ in the 
overall order. This difference — divided by the number of configurations 
involved — is the unfolded level spacing. 

Random Matrix Theory (RMT) predicts three shapes for the statis- 
tical distribution of these level spacings. They refer to the three possible 
patterns of spontaneous chiral flavour symmetry breaking (for Nf flavours) , 

SU{2Nf) S0{2Nf) orthogonal 
SU{Nf) ® SU{Nf) SU{Nf) unitary 

SU{2Nf) Sp{2Nf) symplectic . (4.2) 

At least in 4d Yang-Mills theory with chiral fermions, this set of patterns 
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unfolded level spacing 

Figure 8: The cumulative density of the unfolded level spacing distribution. 
We show the curves corresponding to the RMT prediction for the orthogo- 
nal, the unitary and the symplectic ensemble, and confront them with our 
simulation data (for our lightest fermion mass, m = 0.01 ). We clearly 
observe agreement with the RMT formula for the unitary ensemble. For 
L = 16 there is still a slight deviation for level spacings > 1.5. As we in- 
crease the size to L = 32, even this deviation practically disappears, i.e. 
the agreement becomes very precise. 

is complete {51]0 They correspond to the real, complex and pseudo-real 
fermion representation of the gauge group. (In the real and pseudo-real 
case, fermion and anti-fermion representations are equivalent, hence the 
unbroken symmetry is enlarged to SU{2Nf).) Ref. [56] elaborated the 
corresponding formulae for the unfolded level spacing distributions. For 
lattice QCD (with chiral quarks) the prediction of the unitary ensemble 
has been confirmed but the case of the 2-flavour Schwinger model 

is theoretically less clear, because there is no spontaneous breaking of the 
chiral flavour symmetry. 

Our results (for m = 0.01, as an example) are shown in Figure [81 We 

^"^An overview of the conceivable types of chiral symmetry breaking with an isomor- 
phic representation by non- unitary groups is given in Ref. (55) . 
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L 


m 


(Ai^i,=o) 


(Al, \v\=l) 


(-^1, \v\=2) 


(Al,|z/|=3) 


16 


0.01 


0.1328(6) 


0.175(2) 






16 


0.03 


0.130(2) 


0.173(2) 






16 


0.06 


0.125(2) 


0.173(1) 






16 


0.09 


0.115(2) 


0.172(2) 






16 


0.12 


0.108(2) 


0.166(2) 


0.216(3) 




16 


0.18 


0.108(2) 


0.166(2) 


0.221(4) 




16 


0.24 


0.109(3) 


0.163(2) 


0.215(4) 




20 


0.01 


0.102(2) 


0.127(2) 






24 


0.01 






0.125(4) 


0.148(6) 


28 


0.01 




0.082(3) 




0.120(5) 


32 


0.01 


0.057(3) 


0.076(3) 


0.084(3) 




32 


0.06 


0.030(3) 


0.054(3) 







Table 6: The lowest non-zero eigenvalue of D]^J^p (after the stereographic 
projection Iji4-i\ ) ) for different masses and lattice sizes, in distinct topological 
sectors. 



see very clear agreement with the RMT formula for the unitary ensemble. 
The tiny deviation that we observe for L = 16 is a finite size effect, which 
is manifestly suppressed as we enlarge the lattice size to L = 32. 



4.2 The microscopic Dirac spectrum 

We now focus on the leading non-zero eigenvalue Ai, based on the statistics 
presented in Tables |2] and [31 The results for (Ai) are given in Table |6l We 
mentioned before that chiral RMT has been worked out for the case of 
a finite condensate S(m — )■ 0) (TOlIII], with successful applications in the 
e-regime of QCD. This is not the situation we are dealing with; here S van- 
ishes in the chiral limit, as eq. fll.3p shows. The situation is qualitatively 
similar for fermions interacting through Yang-Mills gauge fields above the 
critical temperature of the chiral phase transition. Also there the under- 
standing of the Dirac spectra is controversial; for numerical studies and 
conjectures we refer to Refs. [581159] . 

In infinite volume, V ^ oo, the chiral condensate is given by the Dirac 
spectrum as 

dX (p : eigenvalue density) . (4.3) 

A + m 



Along with the prediction quoted in Section 1, S oc m^/^, this suggests [19] 

p(A>0)cxAi/S (4.4) 
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in contrast to the Banks-Casher plateau [H] that one obtains in the standard 
setting (with E(m — )■ 0) 7^ 0). In that case, the density for the re-scaled 
small eigenvalues XillV is scale-invariant (at fixed mHV) |60]. In our case, 
the generic relation (Aj) oc [^p(A>0)]~^ suggests that the parameter 

(. = XiV^/^W<^ (for small Xi ; V = L^) (4.5) 

should adopt this role, at fixed /i^ = mV^^^W^ — or simply at small m. 
is a constant of dimension [mass]^/^, which is (in this context) analogous 
to S in the standard setting. 




Figure 9: Cumulative densities of XiV^^^ oc Q, fori = 1...4, at mass 
m = 0.01 and topological charge u = (on the left) resp. |z/| = 1 (on the 
right). We see that Q deviates from scale-invariance. 

Hence we probed the corresponding finite-size scaling, but it is not 
confirmed. This is illustrated in Figure [9] for our lightest fermion mass, 
m = 0.01, in the sectors of topological charge z/ = and |z/| = 1. As a 
quantitative measure, the integrated variance is given — and compared 
to better approaches — in Table [71 Note, however, that the derivation of 
relation (14.41) may be invalidated by inserting an explicitly mass-dependent 
spectral density, p(A,m), in eq. (14.31) . 

Next we consider another scenario, with a reduced exponent of V in 
the re-scaling factor. Now we assume the scale-invariant variable to be 
Zi = XiV'^^^Wz {Wz of dimension [mass]^/'^). This scenario is motivated 
by the fact that it belongs to a theoretically well explored universality class: 
it corresponds to p(A>0) oc A^/^, which is the spectral density obtained 
by RMT in the Gaussian approximation. There is a precise prediction for 
the corresponding spectral density in terms of Airy functions [58] . 

PAiry(^) OC Z[M{-Z)f + [Ai'(-Z)]2 (~ v^/tt at Z > 1) . (4.6) 
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Figure 10: Eigenvalue histograms for Zi oc XiV'^^^ , at m = 0.01, z/ = 
compared to the spectral density pAiry in eq. which RMT predicts 

in the Gaussian approximation. There is no convincing support for this 
scenario. 

Figure [lO] compares the Airy function density (14. 6p to the histograms 
that we obtained in various volumes at m = 0.01 and z/ = 0. Our data 
exhibit a far more marked wiggle structure, hence the agreement is not 
really convincing. 

The finite size scaling quality of Zi is also captured in Table [7l It is 
significantly better than the one of (i, but still not fully satisfactory. 

Let us finally proceed to the most successful approach, which was iden- 
tified empirically. It turns out that our data are in excellent agreement 
with a scale-invariant parameter 

Zi = XiV^^^W;, {W^ : constant of dimension [mass]^/"^) , (4.7) 

which implies a microscopic spectral density p(A) oc A^^^. The scale invari- 
ance of zi . . . is illustrated in Figure [HI again for our lightest fermion 
mass, m = 0.01, in the sectors of topological charge = and = 1, 
which can be compared directly to Figures HI A quantitative confrontation 
with the previous two ansatze is included in Table [71 

We also tested the behaviour if the re-scaled mass is kept approximately 
constant, as an alternative to just keeping m small. In Figure [T^ we com- 
pare {(i) for different lattice sizes, L = 16 and 32, again in the sectors 
1 1^1 = and 1, for //^ ^ const. We add the corresponding test with (zj) and 
fj^z = mV^^^Wz ~ const., which reveals again a clearly superior finite size 
scaling. 

Our data favour this last scenario unambiguously. A hint for a possible 
explanation can be found in Ref. [5], which introduced the dimensionless 




1 2 3 4 5 6 7 
z 
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index 





1 


power p 


3/4 


2/3 


5/8 


3/4 


2/3 


5/8 




1.812 


0.266 


0.026 


2.078 


0.258 


0.066 




0.159 


0.036 


0.007 


0.204 


0.055 


0.019 



Table 7: A measure for the agreement between the cumulative densities of 
the re-scaled eigenvalues Si = XiV^, at m = 0.01 and |z/| = 0, 1. We 
show Vaiint = EJ=i/ dSi[Y,t^{PLf,isi) - Pm(si)) ^(max - 1), the inte- 
grated variance, where the sum over k includes Lk = 16, 20, 32 at u = 
(max = 3), and Lk = 16, 20, 28, 32 at |z/| = 1 (max = 4/ Pm{si) is the 
mean value over the volumes involved. The quantity Var|^^^ is normalised 



by dividing through the relevant interval in s, where 0.01 < Pmis) < 0.99. 
We confirm that the power p = 5/8 yields by far the least variance, i.e. the 
best agreement. 
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Figure 11: Cumulative densities of X-jV^^^ oc Zi, for i = 1 . . . 4, at mass 
m = 0.01 and topological charge 1^ = (on the left) resp. |z/| = 1 (on the 
right). In contrast to Q and Zi, the variable Zi is scale-invariant to an 
impressive precision. 



parameter 

/ := V2mL^/^/{(3'Kf/^ (4.8) 

to distinguish different regimes. The aforementioned behaviour S oc rri^/^ 
is expected for / ^ 1, whereas I <^ 1 <^ ILj^pjf^ imphes S oc vnL. For 
m = 0.01 we are in an intermediate regime, I = 0.5 ... 1.3 (and 2L/ y/iiP = 
8.1 .. . 16.2), which renders our exponent in S oc m^/^ plausible. 

However, a precise explanation for this behaviour remains to be worked 
out. In particular in the framework of RMT — extended to this extraor- 
dinary setting — this might be feasible, but it is far from trivial \61\ . 
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Figure 12: Finite size scaling for {(i) at oc mV'^/'^ ~ const, (left) and 
(zi) at Hz oc mV^f^ ~ const, (right). These plots confirm again that zi 
performs much better as a scale-invariant variable. 

4.3 Eigenvalues in the bulk of the Dirac spectrum 

At last we take a look at Aio as one of the bulk eigenvalues, and we find 
an optimal finite size scaling for Aioi^^'^^, see Figure [13] (left). The plot 
on the right shows that this factor works well also for the re-scaled full 
cumulative density (including all eigenvalues up to the considered value). 
Based on the fact that the spectral cutoff A^ax = 2 is fixed in any volume, 




14 15 16 17 18 19 5 10 15 20 25 



Figure 13: For Aio, a low-lying bulk eigenvalue, the scaling factor is shifted 
to L^'^^. The re- scaled full spectral cumulative densities in different volumes 
(form = 0.01, 1^ = 0) agree well, and turn into the bulk behaviour p(A) oc A 
(resp. Pcumuiative(A) OC X"^ ) , which is expected in d = 2. 

it is now tempting to speculate that the volume factor for a good finite 
size scaling gradually decreases from V^^^ . . .V^. However, considering 
eigenvalues above the regime shown in Figure [131 but below the cutoff 
regime, we could not find any consistent scaling factor. Indeed there is no 
compelling reason for such a factor to exist. 
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5 Topological summation and susceptibility 



The chiral condensate E can be measured using the complete Dirac spec- 
trum, 

where we still refer to the eigenvalues Aj after the projection (14. ip . Here 
we do not need any assumptions from RMT or the e-regime. We want to 
investigate the link to the analytical formula for S(m) in eq. fll.Sp . 

However, since there are only few topological transitions in the HMC 
history (cf. Table [2]), we can only measure 

Sm = -(V'V')Im , (5.2) 

i.e. the chiral condensate in separate topological sectors. Hence an appro- 
priate summation has to be performed. This challenge is generic for HMC 
simulations with dynamical light fermions on fine lattices, so it is relevant 
— in particular in view of QCD — to explore such topological summations. 
Here we encounter an interesting test bed to probe various methods for this 
purpose. 

For very light fermions, the dominant contribution to Sj,^o is due to 
the zero modes. Hence a suitable notation is 

S,, = + £1^1 , eo > ^^1 > > ■ ■ ■ > . (5.3) 
mv 

The inequalities at the end correspond to a general property of stochastic 
Hermitian matrices (such as 73-DovHf): the presence of zero modes pushes 
the small non-zero eigenvalues to higher (absolute) values. 

Our results for the direct measurement of are given in Table El The 
hierarchy anticipated in relation (15. 3 p is consistently confirmed. In addition 
we observe the inequality 

e\,\{Vr) > e\,\{V2) if V^>V2 (5.4) 

to hold. If the volume is enlarged, the non-zero eigenvalues reach out to 
smaller values. We see that this effect supersedes the pre-factor 1/V in 
eq. (15. ip . so that e\y\ increases. The validity of inequalities (15. 3p . (15. 4p is 
illustrated in Figure [TH Moreover, we recognise a smooth mass and volume 
dependence of the terms e\y\. 
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L 


m 


(S.=o) 
(e n) 


(S|.|=i) 

\^\u\ = l/ 


(S|H=2) 

\t-\u\=2/ 


(S|.|=3) 


(^^1=4) 

(£1 1 a) 


16 


0.01 


0.01273(8) 


0.39968(5) 
U.UUyU£)(^£) ) 








16 


0.03 


0.0374(4) 

n nQ7/i (A \ 
U.Uo ( 4(^4 ) 


0.1573(2) 

n no7i (o\ 
U.Uz 1 ^[^) 








16 


0.06 


0.0713(8) 

U.U 1 lOl O I 


0.1174(2) 








16 


0.09 


0.0985(8) 
0985('8) 


0.1182(3) 
0748f3) 








16 


0.12 


0.1174(6) 
0.1174(6) 


0.1275(4) 
0949f4) 


0.1468(4) 
0817('4) 






16 


0.18 


0.1434(3) 
1434('3) 


0.1469(3) 
1252(3) 


0.1548(4) 
0.1114(4) 






16 


0.24 


0.1633(4) 
0.1633(4) 


0.1652(2) 
0.1489(2) 


0.1697(4) 
0.1371(4) 


0.1758(9) 
0.1270(9) 


0.1841(7) 
0.1190(7) 


20 


0.01 


0.0141(3) 
0141f3) 


0.2608(3) 
0108('3) 








24 


0.01 






0.3572(2) 
0.0100(2) 


0.5296(2) 
0.0088(2) 




28 


0.01 




0.1408(2) 
0.0132(2) 




0.3923(2) 
0.0096(2) 




32 


0.01 


0.0181(5) 
0.0181(5) 


0.1112(2) 
0.0135(2) 


0.2073(4) 
0.0120(4) 






32 


0.06 


0.0883(7) 
0.0883(7) 


0.093(1) 
0.077(1) 









Table 8: Results for the directly measured chiral condensate at different 
masses and lattice sizes, in distinct topological sectors. We observe full 
agreement with inequalities Ii5.3\) and ( [5.^| ). As m increases at fixed L, the 
dominant role of the zero mode contributions to S^^q diminished. As L 
increases at fixed m, however, e\y\ is enhanced. 



5.1 Gaussian evaluation of the topological suscepti- 
bility 

In this method we assume a Gaussian distribution of the topological charges, 
so that S can be written as 



exp 



E"=-ooexp(^-2#^j 
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Figure 14: The terms 6\u\ introduced in eq. Ii5.3\) . which represent the contri- 
bution of the non-zero modes to the chiral condensate in a fixed topological 
sector. These plots reveal a smooth and monotonous dependence on the 
fermion mass (on the left, at L = 16) and on the volume (on the right, at 
m = OM). 



where Xt is the topological susceptibihty. At least in QCD the charge 
distribution is indeed Gaussian to a good approximation (see for instance 
the index histograms in Refs. p^). A high statistics study only found a 
tiny deviation, which tends to vanish in the large volume limit [62]. 

If we have data in the sectors up to charge Q, i.e. Q is the maximum 
of the simulated sectors we insert the measured values Eq . . . Eg, with 
the maximum and minimum of the statistical error bar. For higher charges 
we make use of inequality f l5.3p to fix the minimal and maximal values as 

E|i/| min 77 ' E|i/| max El^l rnjn ~l~ Bq . (5.6) 

mV 

In the cases L = 24 and 28 we do not have data for all the sectors with 
1 1^1 < Q. Here we employ in addition inequality (15. 4 p and the results in the 
next smaller (larger) volume to fix S|i,|^nim (X'\u\,ma.x) ■ 

If we insert some susceptibility Xt into eq. (15. 5p . we obtain a value for 
S. Its uncertainty is modest, because most sectors, which have not been 
measured, have exponentially suppressed probabilities p(|z^|). By requiring 
agreement with the theoretical prediction of Ref. [6] (given in eq. (II. 3p ) we 
now determine Xt] the results are given in Table IH] and Figure [TBI 

Table |H] shows that only in the case of light fermions the prediction 
of Ref. [6] can be reproduced. For L = 16 and m > 0.12 all the Si, are 
larger than the prediction, hence there is no way to reproduce it with a 
weighted sum. In the case L = 32, m = 0.06 the value of (S;^=o) is too 
close to the prediction for S to extract a sensible result for Xt ■ At L = 16, 
m = 0.09 a value for Xt which achieves this can still be found, but the result 
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L 


m 


S of Ref. [S] 


Xt 


y) 


16 


0.01 


0.04888 


0.0006586(5) 


0.1686(1) 


16 


0.03 


0.07050 


0.00117(1) 


0.299(3) 


16 


0.06 


0.08883 


0.00159(8) 


0.407(20) 


20 


0.01 


0.04888 


0.000500(2) 


0.200(1) 


24 


0.01 


0.04888 


0.000408(16) 


0.235(9) 


28 


0.01 


0.04888 


0.000367(15) 


0.288(12) 


32 


0.01 


0.04888 


0.000341(4) 


0.349(4) 



Table 9: Results for the topological susceptibility Xt based on the method 
described in Subsection 5.1, which assumes Gaussian charge distribution 
and the chiral condensate according to Ref. 16]. On the L = 16 lattice this 
method does not work for m > 0.09, due to the latter assumption. 




0.005 



Figure 15: The topological susceptibility xt determined by the method de- 
scribed in Subsection 5.1. On the L = 16 lattice we see that Smilga's for- 
mula cannot be valid at m > 0.09, hence also this method is not applicable 
anymore. For the smaller masses our results for the susceptibility suggest a 
behaviour xt oc y/rn. It is compared to the conjecture in eq. ( f5. 8\} . based on 
Refs. [IMU^ . The plot on the right refers to m = 0.01. An infinite volume 
extrapolation, which assumes Xt to be consistent with a linear dependence 
on l/V, leads to Xt = 0.0002263(53). 

does not make sense either, since it is below the Xt values determined at 
m = 0.03 and 0.06 (see Figure ITSl) . We conclude that Smilga's formula 
is not applicable at m > 0.09, which assigns an explicit meaning to his 
assumption m <^ 1/v^ — 0.45. 

Figure [15] (on the left) shows the consistency of our results with the 
behaviour 

Xt{^) oc \fm . (5.7) 
Alternative results (with quenched configurations and re-weighting) were 
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given in Ref. [TH]. On the theoretical side, Ref. [03] conjectured for Nj 
degenerate flavours in the large volume limit 



- = + — , ~ 0.160 g , 5.8 

Xt Sl^^m xt,q 

where S^-^^ := SAr^=i(m = 0) (cf. Section 1), and Xt,q = oo) is the 

quenched value. Actually this mass dependence of Xt was conjectured in 
the framework of QCD. If we apply the same form in the Schwinger model, 
and insert the value Xt,q ~ 0.023(7^ = 0.0046 from Ref. [TH], we obtain 
Xt ~ 0.00033. Figure[T5](on the left) shows a trend towards agreement with 
this conjecture as we increase the mass at L = 16, thus suppressing the 
finite size effects. In particular at m = 0.06 eq. (15.81) implies Xt — 0.0015, 
close to our value in Table M 

Also when we increase the volume at m = 0.01, our results attain the 
vicinity of this prediction: at L = 32 we obtained xt ~ 0.00034, next to 
the predicted value {xt ~ 0.00033). 

If we try an infinite volume extrapolation, however, the plot in Figure 
[15] on the right suggests a linear dependence of Xt on l/V; that assumption 
leads to a smaller value of XtiV = oo) = 0.0002263(53). 



5.2 An approximate formula for the topological sum- 
mation 

The procedure of the previous subsection is robust for light fermions (bar- 
ing finite size effects on S). However, since it uses the analytic result for 
E as an input to determine Xt, it does not evaluate S itself from the nu- 
merical results in distinct topological sectors. In Subsections 5.3 and 5.4 
we will test a method to do so, following Ref. [M]- For convenience we 
re-derive here in a concise form the formula for an approximate topological 
summation that was given in Ref. [6lj for the pion mass; we generalise it to 
arbitrary observables. This consideration follows the lines of Refs. [nHE5|, 
pointing out in particular which assumptions are involved, so the subse- 
quent subsections can refer to them. 

First we assume the fermion field to be integrated out. Thus we refer to 
an effective gauge action S'eff[f/], which keeps track of the fermion determi- 
nant. So the partition function takes the form Z = J T>U exp{— S'eff[f/]}. 
Next we introduce a Kronecker 5 as a filter of gauge configurations \U] with 
a specific topological charge z/, 

5.,u[u] = :r f de exp{ie{u - up])} . (5.9) 
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This formulation involves the vacuum angle ^, and it allows us to write a 
"partition function" restricted to one topological sector as 

Z, = j Pf/e-^^«[^] 5,,,^u] = ^f d9e'''Z{e) , (5.10) 

where Z{9) = J VU exp{— [t/] — i^^z/[[/]} is the (complete) partition 
function for a general vacuum angle (and Z = ^(0)). 

If some observable O is measured only in one topological sector, the 
corresponding expectation value is given bj^ 



a = ^ y I?t/e-^^«[^]5.,.[t;]0[[/] = 1^ f dee'''Z{e)0{e) . (5.11) 
The relation 

-Z"{e)\e=o = {iy') = Vxt (5.12) 

suggests 

Z{9) = Z expi-VxteV2) . (5.13) 
Inserting this ansatz into eq. fl5.10p leads to 

Y = ^f_ d0 ( - + i^'^) ' (5.14) 

where we recognise the stationary phase 

Vxt {v^) 

As an approximation, we extend the bounds in the integral (15.14p to ±oo. 
For this step, it is favourable if Vxt is large. Then we obtain 

^" ^ exp f - . (5.16) 



Z y^2WYt V 2Fx 

This formula is obviously consistent with an integral approximation for 
Z = J2u ^1^' which is again best justified for large Vxt- 

If we insert ansatz fl5.13p into eq. (15. lip , we obtain another integral with 
the stationary phase 6s,u- By repeating its approximation as a Gaussian 
integral, and employing relation fl5.16p . we arrive at 

^ P d9O{e)e^v[-^^{0-9s^^f) . (5.17) 



^"'^We use the notation O, rather than (O), because in the following this will be more 
practical for indicating the dependence on 9. 
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Since we assume only small \9 — ^^.z^l to contribute significantly to this 
integral, we may also approximate 

o{e) ^ (5(0,,,) + ]-o"{e)\e^^{e - • (s-is) 



Let us further assume \9s\ to be small, so we can replace the first term in 
this formula for 0{9) by 

(5(M=^(5(o) + i(5"(0)|oC. (5.19) 

Hence a further property, which is favourable for our approximation, is a 
small topological charge in addition to a large value of Vxt- 

The approximation f l5.19p also implies O"{6)\0_. ~ 0"{6)\q. So we can 
express the (numerically measurable!) restricted expectation value as 



a ^ f d9 [O + \d"[{e - 9,,,f + 0,%]] exp ( - ^{9 - e,,,f) , 

(5.20) 

where O and O" are taken at = 0. Extending once more the boundaries 
to ±cxD leads to the final form 

0.^0 + \0" ^ (l - ^) . (5.21) 
2 Vxt V Vxt) 



This is the same structure as Ref. [M] obtained for the pion mass. It is con- 
sistent that the limit Vxt renders all topological sectors equivalent, 
so that all Oi, coincide with O. 

The numerical measurement with few topological transitions provides 
results for the left-hand-side of eq. f l5.2ip . On the right-hand-side O, O" 
and Xt are unknown. We are interested in O and Xt-, and measurements of 
Oi, in various topological sectors and volumes allow in principle for their 
evaluation (as far as the above approximations make sense). Up to now, 
this intriguing and possibly powerful technique has not been tested with 
simulation data. This will be pioneered in the next two subsections. 

5.3 Topological summation of the chiral condensate 

Now we insert the chiral condensate S as our observable O. As our input 
we have data for some Sj^, i.e. for the left-hand-side of eq. (]5.2ip . but on 
the right-hand side S, S" and Xt are unknown. In view of their evaluation, 
it is convenient to re-write eq. f l5.2ip as 

S" S" A 

A := , B := ^ , Xt = — ■ (5.22) 

2xt 2xr 5 
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Based on data from different topological sectors in a fixed volume V we 
can only evaluate B; it is not possible to obtain E and x without including 
different volumes. 

By considering two sectors with charges = and £ (and k i) at 
fixed V and m, a result for B is obtained as 



1 

V 



B, 



kl 



V 



1 



+ V 



(5.23) 



If more than two values have been measured, the approximate agree- 
ment between the emerging represents a consistency condition. In 
Table dU] we give corresponding results for B^/fV, derived from the data 
in Table E 



m 


L 






B2,i/V 


B^,^/V 


B,,olV 


0.01 


32 


95.33(55) 


48.44(16) 


32.80(15) 






0.12 


16 


2.59(18) 


1.882(50) 


1.647(48) 






0.18 


16 


0.90(11) 


0.730(32) 


0.674(43) 






0.24 


16 


0.49(11) 


0.410(36) 


0.384(38) 


0.356(28) 


0.333(13) 



Table 10: Results for the term Bk//V obtained in each case from two values 
of E|i,| in a fixed volume V and at a fixed mass m, according to eq. (EI 



For small fermion masses, the results for i?^/ are dominated by the semi- 
classical term l/[m(/c -|-£)], and thus strongly dependent on the topological 
sectors involved. This feature is suppressed, however, when m increases; in 
particular for L = 16, m = 0.24 we do observe approximate agreement for 
different choices of k and in striking contrast to the semi-classical result. 
This similarity, which is illustrated in Figure [ini is a remarkable quantum 
effect, due to the fluctuation terms e\y\. Still we observe systematically that 
the Bk^g values decrease if higher topological charges are involved. That 
also reduces the reliability of our assumption of a small \Os,v\ — which 
favours the approximation fl5.19p — so we consider Bi q most reliable. 

In order to proceed, i.e. to get access also to A and thus to x ^'^^ ^i 
we have to confront data from different volumes. Considering two volumes 
Vi and V2, but keeping m and fixed, eq. f l5.22p implies 

^UV.) - EH(Vi) = -4 (^ - ^) + (^ - ^) . (5.24) 

Most convenient is the sector z/ = 0, where we can read off A without 
involving the variable S. At = 1 it is most obvious to insert B as 

..(^0-E.(^.)^^(4-i)+^-^. (5.25) 
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numerical results i — ^ 
semi-classical * 



2 3 
k 



Figure 16: The terms B^q/V for k = 1 . . .4, at m = 0.24, L = 16. For the 
numerical result we observe a remarkable stability in k, in contrast to the 
semi- classical values. 

which is identical to the result obtained from So(Vl), So(^2)- This is 
expected to be the most reliable value for A. Our results obtained in this 
way are listed in Table [TTl 



m 


m = 0.01 


m = 0.06 




(20, 16) 


(32, 16) 


(32, 20) 


(32, 16) 


A 


0.974(221) 


1.833(173) 


2.626(383) 


5.803(363) 



Table 11: Results for the term A in eq. Ii5. 2S\) . obtained from Sq or Ei in 
two volumes, at a fixed mass m, according to eqs. ( [5. 24\ ) and / I5. 25\) . 

The idea is to use these results in both volumes involved. This is sensible 
if A is approximately constant in the volume, but that is not confirmed in 
our data set for m = 0.01. Since the assumptions tend to hold better for 
larger volumes, we use the value of A obtained in (Li,L2) = (32,16) to 
determine S in L = 16, and A from (32, 20) for E in L = 20 and 32. This 
leads to 

S ~ 0.0199(7) (L = 16) , S ~ 0.0207(12) (L = 20 or 32) . (5.26) 

These results for S are nicely consistent, but more than a factor of 2 below 
the value expected in infinite volume. Hence for m = 0.01 this method 
works in an intrinsically consistent way, but the results for S are reduced 
by strong finite size effects. 

Similarly, if we now evaluate Xt by referring to eq. f l5.22p we obtain 
values between 10~^ and 10~^, i.e. well below the results in Subsection 
5.1. Note that the applicability of the method used here — based on the 
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approximations in Subsection 5.2 — is indeed questionable for m = 0.01, 
since we are dealing with a small value of Vxt- 

So the mass m = 0.06 is more promising. We cannot test the volume 
independence of A from our data, but based on L = 16 and 32 we obtain 

m = 0.06 : A = 5.80(36) , S = 0.0940(16) . (5.27) 

Let us focus on the size L = 32 and involve Bi q, which leads to 

Xt = 0.00118(30) . (5.28) 

The method that we used in Subsection 5.1 to evaluate Xt does not work 
in this case as we mentioned before (due to the sizable uncertainty of Si it 
basically just constrains Xt < 0.03). However, if we now insert the preferred 
value Xt — 0.00118 and evaluate S by topological summation, we arrive at 

S = 0.0883(69) . (5.29) 

This result is in excellent agreement with the analytical predictions based 
on Ref. P, S = 0.08883. 

We conclude that some consistency tests are passed well. However, for 
m = 0.01 the ultimate results for xt and the E (summed over all topological 
sectors) seem to be affected by strong finite size effects. 

The situation improves as we proceed to m = 0.06. In this case the 
correlation length (in infinite volume) is only 4.2 lattice spacings, so that 
the finite size effects are suppressed quite well, and we arrive at sensible 
results for xt and S. 

5.4 Topological summation of the "pion" mass 

For practical reasons, it is favourable to consider the decay of a current 
correlation function for measurements of the "pion" mass, rather than the 
pseudoscalar density [7] . In this way we obtained the results in Table [T^ 
Again we can test the topological summation for the fermion masses 
m = 0.01 and 0.06. Here we proceed in a manner different from Subsection 
5.3: if we insert the data of Table [T2] into the formula 

M^,|,| = M^-^ + , (5.30) 

the unknown terms M^-, A and B are over-determined. We choose the 
optimal values for them by a least square fit. 

Let us start again with m = 0.01: if we include all 11 sector that 
we simulated, we obtain Af^ = 0.115(3), and if we skip the 4 sectors with 
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L 


m 










^,jthcory(^ = oo) 


16 


0.01 


0.041(1) 


0.271(4) 






0.071 


16 


0.03 


0.123(5) 


0.275(4) 






0.148 


16 


0.06 


0.214(6) 


0.310(3) 






0.235 


16 


0.09 


0.302(5) 


0.358(2) 






0.308 


16 


0.12 


0.359(5) 


0.413(2) 


0.494(5) 




0.374 


16 


0.18 


0.498(4) 


0.525(2) 


0.589(5) 




0.490 


16 


0.24 


0.631(6) 


0.648(2) 


0.700(3) 




0.593 


20 


0.01 


0.038(2) 


0.209(4) 






0.071 


24 


0.01 






0.257(8) 


0.32(1) 


0.071 


28 


0.01 




0.146(4) 




0.25(1) 


0.071 


32 


0.01 


0.05(1) 


0.160(8) 


0.192(6) 




0.071 


32 


0.06 


0.23(1) 


0.232(7) 






0.235 



Table 12: The "pion" masses measured in various volumes, at different 
fermion masses, in the topological sectors = 0...3. The last column 
displays the (infinite volume) prediction of Ref. (cf. eq. (T^). 



|z/| > 1 (which are problematic for our approximations) the result is slightly 
reduced to M^^ = 0.113(4). Thus we confirm the observation of Subsection 
5.3 that the summation formula runs into trouble for this light mass. Ac- 
tually this is not surprising: we recall that the correlation length for this 
fermion mass (in infinite volume) is ^ ~ 14, cf . Figure HJ hence strong finite 
size effects are expected, and they generically enhance the mass gap. 

Therefore, we now skip the smallest volumes. We need at least two 
volumes to determine the term A, so we restrict the consideration to L = 28 
and 32, and we include < 2; thus we are left with four sectors. This 
small number of input data causes a large error, but the optimal value of 
the least square fit moves very close to the theoretical prediction (see Table 

m = 0.01 : = 0.073(25) . 

For m = 0.06 we only have four sectors to deal with, but the finite size 
effects are much less severe, and we arrive again at a sensible result, 

m = 0.06 : = 0.232(8) , Xt = 0.0007(4) . (5.31) 

As we already saw in Subsection 5.3, this method is not very useful for the 
determination of Xu it is always plagued by a large uncertainty, as in the 
example given here (for m = 0.01 this is even worse). The three results 
that we obtained for at m = 0.06 (in Subsection 5.1, 5.3 and 5.4) differ 
within the same magnitude^ 

^^The result in Table |9] is larger, but only based on L = 16. 
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Regarding the "pion" mass, we confirm that this method can work, if 
it is apphed in an appropriate regime. Hence this approach may in fact be 
promising for further apphcations, including QCD, if finite size effects are 
under control — in particular the term (z/^) should not be too small. 

5.5 Correlation of the topological charge density 

Members of the JLQCD Collaboration proposed an alternative method to 
evaluate Xt even in one single topological sector based on the topological 
charge density, pt{x). Ref. |65] derived the "model independent formula" 

l.lZo {Pt{x)ptm\H^-f + ^,{iy' + ^J+0{V-') , (5.32) 

which captures even a possible deviation from the Gaussian charge distri- 
bution by a non-vanishing kurtosis C4 = {{i^^) — 3{i''^))/V However, this 
term is known to be tiny, so we neglect it in the following. 

The issue is to search for a plateau value of the charge density cor- 
relation at large distances, which differs from the constant i/^/V"^ (in the 
sector with topological charge This shift A = —Xt/V is negative 

because fluctuations in pt{x) have to compensate. However, A tends to 
be small and hard to resolve numerically. For the absolute value |A| it 
is favourable if m increases (although we are actually interested in quasi- 
chiral fermions), but it is unfavourable if V increases (although we need 
access to the asymptotic value in eq. f l5.32p . and 0(y~^) should be negli- 
gible). Moreover, in analogy to the approximations that lead to eq. (15.211) . 
the derivation of eq. fl5.32p involves the assumptions 

(z/^) = Vxt is large , is small . (5.33) 

Our estimates for xt suggest that we do not have any setting which is 
really adequate for the flrst of these two conditions. Sometimes, however, 
such approximations apply reasonably well even if the assumptions do not 
strictly hold (this is the experience in QCD with simulation data matched 
to formulae of the e- or 5-regime, and in the preceding two subsections to 
some extent). 

Since our conflgurations are smooth, it is not problematic to use the 
naive lattice formulation of the topological charge density, pt = ei2-Fi2- 
(There is a cleaner formulation for Ginsparg- Wilson fermions [25], but it 
is tedious in practice.) As examples, the results for L = 16 in the topolog- 
ically neutral sector for three masses are shown in Figure [T71 In particular 

-"^^A variant of this method, which uses the Ty'-correlation of the pseudo-scalar density, 
has been apphed to two-flavour QCD in Ref. i66j . 
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Figure 17: The topological charge correlation function for L = 16, u = 
and m = 0.01, 0.03 and 0.06. 

we see that the correlation over short (but non- vanishing) distances is neg- 
ative, since a given hnk variable contributes with opposite signs to pt in 
its adjacent plaquettes. The jackknife errors for the quantity {pt{x)pt{0))\o 
are typically in the range (1 ... 2) ■ 10~^. 

In the most promising case, m = 0.06, we evaluated 741 topologically 
neutral configurations and we have identified the magnitude for Xt ~ 0.0015 
which implies A ^ —6 ■ 10"^. In order to safely resolve this shift we would 
therefore need about 50 000 to 100 000 configurations. We conclude that 
this method requires a very large statistics, so we cannot apply it. 



A 



0.0001 



-0.0001 



-0.0002 



6 Conclusions 

We presented a study of the 2-fiavour Schwinger model with dynamical 
chiral fermions. We applied the overlap Hypercube Fermion (HF), and 
confirmed its features as a promising formulation of a chiral lattice fermion. 
This is manifest by its excellent locality and scaling behaviour. Moreover, 
our study allowed us to explore a number of conceptually interesting and 
relevant issues, and to test new methods to handle them. 

This is certainly of interest, even in a toy model, since simulations with 
dynamical overlap fermions have been explored only poorly so far. 

For our simulation we designed a specifically suitable variant of the 
Hybrid Monte Carlo (HMC) algorithm, demonstrated its correctness and 
tested its efficiency. It is sufficient to insert a low polynomial approxima- 
tion for the overlap operator to compute the force term, while the high 
precision overlap operator is employed in the accept/reject step. It is an 
open question if this strategy — perhaps with further refinements — can 
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be carried over to the simulation of QCD with dynamical chiral quarks. 

Next we discussed the spectrum of the Dirac operator. Random Matrix 
Theory is not yet worked out for this type of model, with a vanishing chiral 
condensate at zero fermion mass. Nevertheless the unfolded level spacing 
density follows the standard RMT formula for the unitary ensemble. 

The prediction E(m) oc rn^l^ [IHS] suggests a microscopic spectral den- 
sity p(A>0) oc A^/^, and the scale-invariant variable C oc XV^^^ (if we 
assume no explicit mass dependence of p). This conjecture does, however, 
not agree with our data. An alternative scenario with p(A > 0) oc A^/^ has 
a theoretical background as well (Gaussian distribution), but the data do 
not really support it either. Instead they favour z oc XV^^^ as the scale- 
invariant variable, and therefore p(A > 0) oc A^''^. 

A hint for an interpretation of this surprising result can be found in 
Ref. [5], which derived the behaviours S oc m}/'^ and S oc mL in two limit- 
ing cases. Regarding the parameter that characterises these extreme cases, 
our settings are in an intermediate regime, which appears compatible with 
the relation E oc rn?/^ that we observed. For a precise theoretical clarifica- 
tion, we hope for the corresponding RMT formulae to be elaborated. 

Direct measurements of the chiral condensate and of the mass of the iso- 
triplet ( "pion" ) could only be performed in fixed topological sectors, since 
the HMC histories contain only few topological transitions. This limitation 
is generic for the simulation of light fermions close to the continuum limit. 
It is therefore a major challenge to explore techniques for the topological 
summation of such measurements. Assuming a Gaussian distribution of 
the topological charges, the data can be used to evaluate the topological 
susceptibility. In order to determine the actual observable, we explored 
approximate summation techniques, which led to sensible results in some 
parameter window, where the term {y"^) = Vxt is not too small. 

Topological summations could become relevant for future QCD simula- 
tions. Nowadays they are carried out with very light dynamical quarks — 
such that the pion mass is close to its physical value — on finer and finer 
lattices. In particular in applications of dynamical overlap fermions |T7] 
the HMC history tends to be trapped in the topologically trivial sector for 
the entire simulation, which endangers ergodicity and does not provide the 
physical result (unless the volume is very large). For Wilson-type lattice 
fermions (which break the chiral symmetry explicitly), the problem is less 
obvious on the currently used lattice spacings, but on still finer lattices 
(with a < 0.05 fm) it is expected to show up as well [22]. Hence topologi- 
cal summation methods are of interest for non-perturbative studies of low 
energy nuclear physics based on first principles of QCD. Our results tell us 
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to be cautious with such summations, but at the same time they provide 
hope for their feasibihty. 
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